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Abstract 

A new model of radiative neutrino masses generated via two-loop diagrams is 
proposed involving a charge 2/3 vector-like quark and a doublet of leptoquark scalars. 
This model predicts one of the neutrinos to be massless and admits both the normal 
and inverted neutrino mass hierarchies with correlated predictions for li — >• £j + 7 
branching ratios. New contributions to CP violation in B s — B s mixing arise in the 
model through leptoquark box diagrams, which can explain the anomalous dimuon 
events reported by the D0 collaboration. These leptoquarks, with masses below 500 
GeV, also provide a natural resolution to the apparent discrepancy in the measured 
values of the CP violation parameters sin 2/3 and €k ■ 
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1 Introduction 



Neutrinos must have tiny masses, so that different flavors can oscillate among one another, 
as observed in experiments. An elegant and natural way to generate the tiny masses is 
through the dimension-five lepton number violating operator C = 0\jM where [f] 

0i = UUH k H l e ik e jl . (1) 

Here L is the lepton doublet and H the Higgs doublet, with i,j = 1,2 being SU(2) L 
indices. The suppression by an inverse power of M, which can be much greater than the 
weak scale, explains the smallness of neutrino mass, which is given by m v ~ v 2 /M, with 
(H°) = v ~ 174 GeV being the Higgs boson vacuum expectation value (VEV). Operator 
0\ is naturally realized through the seesaw mechanism wherein right-handed neutrinos, 
which are singlets of the standard model (SM) gauge group with large Majorana masses, 
are integrated out [2]. The effective mass scale M should be of order 10 14 GeV in order to 
generate neutrino masses of order 0.1 eV, as indicated by neutrino oscillation experiments. 
Such a large scale of M would however make this mechanism difficult to test directly in 
experiments such as the ones pursued at the Large Hadron Collider. 

An alternative method for inducing naturally small neutrino masses is the radiative 
mass generation mechanism [3-7]. This scheme posits that the dimension 5 operator 0\ of 
Eq. (1) is absent, or is highly suppressed, so that neutrino masses remain zero at the tree 
level. Lepton number violation arises through effective operators with dimension d > 5, 
typically containing charged fermions as well as the neutrino fields. These operators can be 
converted to neutrino mass, but only through loop diagrams, wherein all charged fermions 
are annihilated. The induced neutrino masses are naturally small, even when new particles 
needed to generate the d > 5 lepton number violating operators have masses in the TeV 
range, owing to chirality and loop suppression factors. 

The simplest set of operators carrying two units of lepton number appropriate for small 
Majorana neutrino mass generation, in the absence of 0\ of Eq. (1), is of dimension seven. 
There are six such d — 7 operators [8] : 

2 = VVL k e c H l e ZJ e kl 

3 = {L i UQ k d c H l e ij e kU UUQ k d c H l e^} 

4 = {L i VQ i u c H k e jk , L i VQ k u°H k e ij } 

O s = ISe c u c d c H j eij . (2) 

Here the generation and color indices have been suppressed. Q, L denote left-handed quark 
and lepton doublets, while u c , d c , e c denote left-handed anti-quark and anti-lepton singlets 
of the standard model. A full list of AL = 2 effective operators through d — 11 is given 
in Ref. [8]. Among the d = 7 operators of Eq. (2), 2 is perhaps the simplest, which 
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can be induced when the scalar spectrum of the standard model is extended to include a 
second Higgs boson doublet and a charged singlet scalar field h ± . This is the well-studied 
Zee model of neutrino masses [3]. In its simplest version, with natural flavor conservation 
in the Higgs sector, this model predicts vanishing diagonal elements of the neutrino mass 
matrix [3,9], which is now excluded by neutrino oscillation data [10]. 

A second widely studied model of radiative neutrino masse generation [6,7] has a purely 
leptonic effective d — 9 operator, Og = L l L? L k e c L l e c e^-e^, suppressed by M -5 . Here 
neutrino masses are induced via two-loop diagrams. This operator can be obtained when 
the standard model is extended to include a singly charged (h + ) scalar and a doubly charged 
{k ++ ) scalar. The resulting model fits the neutrino oscillation data well, and also predicts a 
host of leptonic flavor violation processes, some of which within reach of ongoing and next 
generation experiments [11]. Operator 0% of Eq. (2) is best induced by scalar leptoquarks, 
as recently shown by us in Ref. [12]. This model leads to consistent neutrino phenomenology 
and interesting flavor effects in both the quark and the lepton sectors [12]. For discussions 
of models based on other operators, see Ref. [8, 13, 14]. 

Operator O3 of Eq. (2) is the main focus of this paper. It has two different SU(2)l 
contractions possible, as shown in Eq. (2). These operators arise in supersymmetric models 
with i?-parity violation. The superpotential couplings W D XLLe c + \'QLd c would gener- 
ate O3 once the SUSY particles are integrated out [5,15]. The QLd c term would induce the 
second contraction of O3 in Eq. (2), the LLe c term would induce O2, while the product of 
QLd c and LLe c would induce the first contraction of O3 . There is an important difference 
in the first and second SU(2)l contractions of O3. In the second contraction, neutrino 
masses are induced at the one loop, while in the first contraction, they arise only at the two 
loop level. (In the second contraction, two neutrino fields appear, while the first has one 
neutrino field and a charged lepton field, which must be annihilated to convert this operator 
to neutrino mass.) The focus of this paper is models which induce the first contraction of 
O3, without inducing other operators that lead to one loop neutrino masses. SUSY with 
i?-parity violation does not fit this requirement, as O2 and/or the second contraction of 
O3 are also induced there. The simplest possibility we have found is to add a vector-like 
charge 2/3 iso-singlet quark to the SM, along with a doublet of leptoquark scalars. The 
induced neutrino mass is of the form 

fghXb v 2 

where f,g,h are dimensionless Yukawa couplings, A& = rrib/v is the 6-quark Yukawa 
coupling, and M stands for an effective mass of the vector-like quark/leptoquark. For 
/ ~ g ~ h ~ 10~ 2 , the mass scale M should be of order TeV, in order to generate m u ~ 0.1 
eV. It is, however, evident from Eq. (3) that M can be as large as about 10 8 GeV, when 
f,g,h are of order one. There are several reasons for considering low values of M, first 
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and foremost being direct tests of the vector quark and the leptoquarks at the LHC There 
are hints of new physics in the B meson system, which can be explained by the new lep- 
toquark scalars and/or the vector-like quark of the present model. The D0 collaboration 
has reported an excess in the same sign di-muon asymmetry in B decays [16], which may 
be a hint for new CP violation in B s — B s mixing. There has also been a tension in the 
determinations of the CP asymmetry parameters sin 2/3 in B meson decay and ex in Kaon 
decay, which may need new physics [17]. The present model, with leptoquark masses below 
a TeV, can explain these anomalies. Furthermore, when this model is eventually embedded 
in a supersymmetric framework, M of Eq. (3) will have to be close to the SUSY breaking 
scale, owing to the SUSY non-renormalization theorem, with the consequence that all loop 
diagrams that generate neutrino masses cancel in the supersymmetric limit. 

This rest of the paper is organized as follows. In Sec. 2 we present the model leading 
to the two-loop neutrino mass generation via 0%, the first contraction of Eq. (2). In Sec. 3 
we obtain the experimental constraints on the model parameters arising from rare process 
in the quark as well as lepton sectors. Here we show how the proposed model explains the 
discrepancy observed by D0 in the CP asymmetry of the B s system. New contributions to 
the CP violating decay B^ — > J/^/Ks are shown to be of the right magnitude to explain the 
apparent tension between sin 2/3 and ex determination. In Sec. 3 we also evaluate the rate 
for neutrinoless double beta decay induced via the vector-scalar exchange mechanism [18]. 
Finally, we give our conclusions in Sec. 4. 

2 Radiative neutrino mass model with vector— like quark 

We wish to generate the operator (L • L){Q ■ H)d c in a renormalizable theory. Here and 
in discussions that follow we use a compact dot product notation for SU(2)l contraction: 
L ■ L = L' l L j eij, Q ■ H = Q^-H^e^, etc. The simplest way to generate this operator, without 
inducing other operators that generate neutrino masses at one loop, is by integrating out a 
charge 2/3 iso-singlet vector-like quark, and a doublet of scalar leptoquarks. These fields 
transform under SU(3) C x SU(2) L x U(1) Y as 



Fermions : £7(3, 1, 2/3) + U%3*, 1, -2/3), Scalars : fi(3, 2, 1/6) = _ 1/3 . (4) 



These particles will have new Yukawa interactions with the SM fermions as well as gauge 
invariant masses given by 



where Q = ir 2 Q*. Here the dots indicates SU{2)i J contraction, as mentioned earlier, and we 
use indices i,j to denote generations. Possible mass terms miU^U, not shown in Eq. (5), 





(5) 
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can be rotated away by field redefinitions. The simultaneous presence of the interaction 
terms gij, hi, fi would lead to lepton number violation by two units, a necessary condition 
for neutrino mass generation. 

We should also specify the scalar interactions that couple the leptoquark Q with the SM 
Higgs doublet H. There is a single non-trivial quartic coupling between these two fields: 

^Zt = M^-H\ 2 (6) 

When the neutral component of the SM Higgs doublet H° acquires a VEV, this quartic 
coupling will generate a mass splitting between u 2 ^ 3 and w _1//3 leptoquarks: 

M 2 ^ 1/3 =Ml M 2 2/3 = Ml — Ml — Xv 2 , (7) 

where v = \[2m w j g ~ 174 GeV. 

The mass matrix for the charge 2/3 quarks, including U, U c fields, that follows from Eq. 
(5) has the form 

M,= ( Y f! ,°A (8) 



fv M 

where (u°, U c ) multiply on the left and (ui, U) multiply on the right. Here Y u is a 3 x 3 
Yukawa coupling matrix, / is a 1 x 3 row vector, and stands for the 3x1 null column 
vector. This mass matrix can be diagonalized by a biunitary transformation 

M u d = UM U V ] (9) 

where U, V are 4x4 unitary matrices. Without loss of generality we choose a basis where 
the 3x3 matrices for the down quarks and charged leptons are diagonal. Thus, the CKM 
matrix will be the 4x3 sub-matrix of the 4x4 matrix V. The charged current interactions 
of the quarks, therefore, become 

Clr = ^u a V ai r(l -K)diW+ + h.c. : (10) 
C%L = ^^n a [(M^) a V m (l- l5 )-V ai (M d ) i a + l 5 )}d t H + + h.c. (11) 

The Greek indices a, — 1 — 4 label generations in the up-quark sector {u\ — u, U2 = 
c, W3 — t, «4 = t'), while the Latin indices i,j — 1 — 3 label generations in the down-quark 
and lepton sectors. Introduction of vector-like quarks U, U c to the SM spectrum will induce 
flavor changing neutral currents (FCNC) in the charge 2/3 quark sector, which are given 
by 



CeX = j-^^[^7 M (l-|sin 2 ^-75)-^A; 4 7 M (l-75)]^ (12) 

4 COS t>w 

CaL = 7r4— u a (M u ) a V aj V^l- l5 )u^H° + h.c. (13) 



2y/2m 



w 
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These interactions can generate tree-level D — D mixing, as discussed in the next section, 
which will strongly constrain the product | V14V24] . 

The 4x4 unitary matrix V can be parameterized as [19] 



V = 
( 



C12C13CU 



-C23C24S12 

-Cl2Ci 3 Si4S24e 



Cl2C24Sl3S23e 

t(l5l4 — $24.) 



i$13 



C34S12S23 - Ci 2 C23C34Si 3 e 
-Cl2Ci3C24Si4S34e l<514 
+ C 2 3Sl2S24S 3 4e 1 ' 524 
+ Cl2S23S24S 3 4Sl 3 e l(,5l3+,524) 



-c 12 c 13 C24C 3 4Sue tSl4 

+Cl2C23S 13 S 3 4e lSl3 
-C23C34Sl 2 S24e l<524 - S12S23S34 
+ Cl2C 3 4Sl3S23S24e 4 (' 5l3+ ' 524 ) 



C13C14S12 



C12C23C24 

-Cl3Sl2Sl4S24e 



C24Sl2Sl3S23e 

(#14 — 624) 



i<5i3 



i$i3 _ 



C12C34S23 - C 2 3C34Si 2 Sl3e 



S<5l3 



— Cl2C23 s 24S34e 



-J&2 



-Cl3C 2 4Sl2Sl4S34e 
+ Si 2 S 2 3S24S34Sl3e i(<5l3+<52 



"- 1 +C12S23S34 
-Ci3C 2 4C34Sl2Sue lSl4 



-Cl2C23C34S24e 



+ C23Sl2Sl3S 3 4e l<513 
+ C 3 4Sl2Sl 3 S23S24e 4 ( ,5l3+ ' 524 ) 



ci 4 si 3 e 



-iS 1 



C13C24S23 Ci 4 S 2 4e 

(<5l3+<524— <5l4) 



-Sl 3 Si4S 2 4e 

C13C23C34 



-Ci3S23«24S34e ' 



-C24Sl3Sl4S34e 



-C13C23S34 
-Cl3C34S23«24e 



-C24C 3 4Si3Si 4 e 



i(5i4-<5i3) 



'id 2 



,(5i4.-S 13 ) 



'524 



C14C24S34 



C14C24C34 



(14) 



where s a/ g = sm6 a p, c a p = cos6 a p. The CKM mixing matrix elements V a i are the elements 
of the 4x3 sub-matrix of V. In terms of the fermion mass eigenstates, Eq. (5) can be 
written as 

JOT = dj {g T ) 3l (^ur 1 / 3 - e^ 3 ) 

- {uJC t uj-^ + iJCW/ 3 ) K V:^^u a + h.c. (15) 
which will be used in our calculations. 

2.1 Two- loop neutrino masses 

By combining the interactions given in Eqs. (6), (10), (11) and (15), one can construct 
diagrams generating neutrino masses. These diagrams arise at the two loop level, and are 
shown in Fig. 1. We have done the evaluation of these diagrams in general gauge, so the 
unphysical Goldstone mode H + also appear in this set. A non-trivial check of the calculation 
is the gauge independence of the induced neutrino mass, which we will show explicitly. Since 
the external neutrinos are Majorana particles, there is another set of diagrams identical to 
the ones in Fig. 1, but with all internal particles replaced by their charge conjugates. The 
sum of these diagrams would make the neutrino mass matrix symmetric in flavor space. 

The induced neutrino mass matrix is proportional to the down quark mass matrix, since 
these diagrams make use of the SM charged currents, which require a chirality flip for the 
d c fields. This is explicitly shown in Fig. 1. The neutrino mass matrix, therefore, can be 
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(H°) 
t 

, — i — 

, ,-1/3 ' y , ,2/3 

t I \ 

-» 1 < X > 1 < )< > 1 — <- 

Vi d% d k u% u a Vj 
(l) 



-1/3/ 



— r - ^ 

' / ,2/3 



(2) 



"1/3, 



-> 1 < X > 



X , ,2/3 
I 

-> — 1 — <- 



(3) 



-1/3 



-1/3 



-1/3 



-) — 1 <)<)>> 




(4) 



-) — 1 < X > i < X > 1 < X > i < 



Vi d% d k \ u c a u a &j e c 3 i Vj Vi d% \ u a 



-» 1 <X) | < 



H+ 



(5) 



(6) 



,,2/3 



.2/3 



J 2 /3 



< x ) j-M — «- 



> I < X > 1 < )<> i < X > 1 — <- 

Vi \ e- ej <f| d k i u c a u a Vj 
\ i 
\ / 

* - ■*- - ' 



> i <x> ' < , — > — 1 — <- 

Vi ve| ej <i| / «i a u i 
\ i 
\ ✓ 

- ' 

H+ 



(7) (8) (9) 

Figure 1: Two- loop diagrams leading to finite neutrino masses in general gauge. 
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:9 2 m b 



hi(V^) 4:a V ak (D d ) k (g T ) k jI aki j + gi k (D d ) k (V T ) ka (V*) a 4 : hjI ak j i , (16) 



written as 

where 3 is a color factor and D d is the normalized down quark mass matrix, 

D d = diag 



m d m s 
m 6 ' m b ' 



The function I aki j is a sum of loop integrals defined as 

3 6 9 

f \^ f(n) , \^ f(n) , \^ f(n) 

1 akij Z_^ 1 ak ~ / ^ 1 aki / j 1 aki ' 



n=l 



n=4 



n=7 



where the integral J*"-* is given by 1 



f(i) , f(2) = 



i 2 w J 7 (2tt)4 (2^)4 fc2 ^ 



/c 2 — £m 2 



(17) 



(18) 



q 2 -M 2 q 2 -m 2 dk (k + q) 2 — Mf (fc + g) 2 - m 2 a ' 



(19) 



f(3) 

1 ak — 



I 



d A k d A q fl+i 



(2tt) 4 (2tt) 4 fc 2 -m 2 



frk-(k + 2q) f i , k 2 - m 2 w 



k 2 



• k 2 - £m 2 



w 



q 2 -M 2 q 2 -m 2 k (k + q) 2 - M 2 (k + q) 2 



(20) 



f(4) 
aki 



d 4 k d A q 



(2^) 4 (2tt) 4 
1 



1 



g 2 -M 2 g 2 -m 2 (fc + g ) 2 -m 2 Q ' 



1 



k 2 — m^y k 2 — m 2 . 



(21) 



r(5) , f(6) _ _ / m e, \ /" d fc A (ft A)ft / 

+ aki ~ \m w ) 7 (2tt) 4 (2^) 4 fc 2 v ^ 



aki ~ aki 



1 



(2tt) 4 (2^) 4 k 2 \ s k 2 -£m 2 w J k 2 



rat 



q 2 -M 2 q 2 -m 2 (k + q) 2 m 2 ' 



(22) 



9(7) 

akj 



-I 



d 4 k d 4 q (fi+ffl 



(2tt) 4 (2tt) 4 k 2 ~m 2 w 



3-C 



k 2 - 



w 



k 2 — ^rrih, J k 2 — m 2 A 



q 2 — m 2 



dk 



(k + q) 2 -m 2 Ua (k + q) 2 -M 2 



2 ' 



(23) 



f(8) ,f(9) 
akj ' akj 



in, 



mw 



J (2^) 4 (2tt) 4 k 2 \ ? 



k 2 - im 2 w 



k 2 



q 2 -m 2 k (k + q) 2 -m 2 Ua (k + q) 2 - M 2 



2 ■ 



(24) 



1 Owing to the unitarity of V, only terms containing m Ua are relevant in generating neutrino mass (see 
Eq. (16)). All terms that are independent of m Ua will add up to zero, and therefore, such terms are not 
written explicitly. 
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It is straightforward to show that all terms containing the gauge parameter £ in Eqs. 
(19)-(24) add up to zero. This means that the neutrino mass matrix elements, which are all 
physical, are gauge independent. An interplay of all diagrams of Fig. 1 is required to see this 
gauge independence, although this can be inferred before doing the momentum integrals. 
Note that the contributions to these integrals proportional to charged lepton masses are 
strongly suppressed as can be seen from Eqs. (21)-(24). Thus, it is a good approximation 
to work in the limit m ei ~ 0. In this limit, the neutrino mass matrix is reduced to a 
rank two matrix with a suppressed determinant det (M v ) <C (0.01 eV) 3 . Thus, we have a 
prediction that the lightest neutrino is essentially massless. For the purpose of evaluating 
these integrals we can also set the down quark masses to zero. Thus, the neutrino mass 
matrix of Eq. (16) can be written as 



(M„) . ~ -g 2 m b [hi{V Y ) ia V ak (D d ) k (g T ) kj + g ik (D d ) k V kc y* 4 h 3 ] I a , 



(25) 



with I n 



Lakij{1Tl dk ~ m e 



~ m e . ~ 



0). The asymptotic behavior of integral I a (in the 



limit Mi = M%) is given by 



-6 In 



Mi 



( ™ Ua \ 2 _ tP_ 
\ Mi J 2 



+ 5 

~ 2 



7r 2 + ^, for Mi > m Ua ,m w , 



ln ' ( ! lif) 2 + T + 6 ' form Ua »Mi, 



(26) 



6-x 



Mi 



\ + 9, for m Ua ~ Mi > m W - 



These expressions are very helpful, especially for analytic approximations of the integrals 
where the internal quarks are light quarks, and also as cross checks of the exact numerical 
calculations. 

The neutrino mass matrix can now be written down: 

/ 



M„ 



m 



1 hi 



1 h 2 



2 h\y 2/11 

V 2 h >. ~ 2 hi 



X 



2h 2 y 1 2/ii 

y 

2/13 ^ 2h 2 y 



lh 
2 h 3 



Hi. _|_ \ 
ho "r 2 /ll \ 



2 h 3 2 /i 2 y 



(27) 



where 



m = 



3g 2 m b 
(16tt 2 ) 2 



hiF\ 

hF 3 



y = 



h 2 F 2 
h 3 F 3 



(28) 



with repeated indices assumed to be summed. By using the unitarity of the mixing matrix 
V, and the fact that m u , m c <^ m t , m t > we have 



v: 4 v ak i a ~ v 3 * 4 v 3k (/ 3 - h) + v; 4 v 4k (h - h] 



(29) 



9 



0.5 1.0 1.5 2.0 2.5 3.0 

My (TeV) 



Figure 2: Plots of the integral functions ^3 — l\ and I4 — I\ versus the leptoquark mass Mi. 
The mass of the vector-like quark is taken here to be 600 GeV, with M\ — M 2 = 60 GeV. 

Plots of I3 — Ii and I4 — Ii as function of the leptoquark mass Mi are shown in Fig. 2, for 
a fixed value of the vector-quark mass of 600 GeV. 

To illustrate the range of parameters allowed by the neutrino mass, let us assume <C 
gj 3 , k = 1, 2, so that only gj 3 contribute to the neutrino mass matrix. If we further assume 
that only the the top quark (among u, c, t) mixes significantly with the vector-like quark, 
i.e., fx, f 2 < then V34V33 ~ -V44V43 ~ hv/M. Therefore, we can write 

Fj * 9j jf (h - h) • (30) 

For normal neutrino mass hierarchy, m — 0.03 eV is needed, which in turn requires ^3^3 — 
10 . This means that for order one values of the Yukawa couplings h^, f$, g^, the mass of 
the vector-like quark and/or the leptoquarks is of order 10 8 GeV. Conversely, if both the LQ 
and vector-like quark have masses of order TeV, and if h% ~ 1, one must have (733/3 ~ 10" 5 . 
In both regimes, lepton flavor violation processes do not strongly constrain the model 
parameters. Interesting new effects will arise, however, if the vector-like quark/leptoquark 
masses are near a TeV, and if some of the Yukawa couplings lie in the range 10 -2 — 1, as 
will be discussed in the next section. 

Although the model predicts the lightest neutrino to be essentially massless, owing to 
the highly suppressed determinant of M u , Eq. (28) does admit both the normal hierarchy 
(NH) and the inverted hierarchy (IH) of neutrino masses. Since the off-diagonal elements 
of M v are uniquely related to the diagonal elements, one can determine the values of hi/hj 
for i < j as 



1 ± 



{M v )jAM v 



(31) 



10 



where (M v )ij are obtained from 

M v = Up MNS (M u ) diSLg Ul MNS . (32) 

Here ?7pmns is the leptonic mixing matrix parameterized as in Ref. [20], while (M^)diag is 
given by 

(M,) diag = diag (0,m 2 e ia ,m 3 ) , for NH, 

(M,) diag = diag ( mi ,m 2 e ia ,0) , for IH. (33) 

Take for example the ratio h\/h 3 . Its value can be determined from Eq. (31), but this 
must match the product {h\/h 2 ) ■ (h 2 /h 3 ), also obtained from the same equation. Now, by 
using the central values of the current neutrino oscillation data, Am 2 ol = 7.59 x 1CT 5 eV 2 , 
Am 2 tm = 2.3 x 1(T 3 eV 2 , sin 2 6 12 = 0.304, sin 2 6 23 = 0.5, and the upper limit on 13 , 
sin 2 < 0.04 [21], one can find the allowed values of h\/h 3 . This result is plotted in Fig. 
3 versus sin 2 13 , both for NH (upper left panel) and for IH (upper right panel). From this 
figure we see that the ratio of h 2 /h 3 has to be of order one for normal hierarchy, while it 
can range from 0.5 to 1000 for inverted hierarchy. In both cases, the value of #13 is allowed 
to range from zero up to the current upper limit. Recently, the T2K experiment [22] has 
reported an indication of nonzero $13, with the best fit value (assuming sin 2 29 23 = 1 and 
5 = 0) being sin 2 2^13 = 0.11 (0.14) for normal (inverted) hierarchy. MINOS experiment 
also finds supporting evidence, although with less significance [23]. The present model can 
accommodate these indications for a sizable 6 13 . 

The couplings hi will mediate £i — > Iff decays (see the next section for detailed discus- 
sions). One has for the ratio of branching ratios, 



BR (/I ->■ erf) 



2 



x BR(r ->■ eV e v T ), (34) 



BR (r ->■ ej) 

where BR(r — > eV e u T ) ~ 0.18 [20]. Eq. (34) has an interesting consequence. As explained 
above, in the NH case, \h 2 /hz\ ~ 1. This means the branching ratio of r — > cannot exceed 
5.3 x 10~ n because of the limit on the branching ratio Br(/i ->■ ej) < 2.4 x 10~ 12 [25]. A 
measurement of BR (r — > ej) near the current experimental limit of ~ 10~ 8 would rule 
out the NH scenario. Of course, for these decays to have significant branching ratios, the 
leptoquarks must have masses not much above a TeV. In Fig. 3, we also show the ratio 
\hi/h 2 \ as a function of sin 2 9 13 allowed in the model for the NH case (lower left panel) and 
IH case (lower right panel). The ratio Br(r — > e7)/Br(r — > /ry) = \hi/h 2 \ 2 in our model, 
which can server as a further test. 



3 Experimental constraints 

The new interactions shown in Eqs. (5) and (15) can induce lepton flavor violation processes 
such as ji — > e7 and /i — > 3e decays. In this section we analyze various such processes and 
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Figure 3: The plot of the allowed value of hi/hj as a function of sin 2 6*13. 



derive limits on model parameters. LHC experiments have set lower limits on the leptoquark 
mass: M\p > 376 (319) GeV for the first generation leptoquarks and M\p > 422 (362) GeV 
for the second generation leptoquarks, assuming branching ratio of 1(0.5) [24]. Our fit to 
neutrino mass suggests that the branching ratio of the leptoquark to muons is about 0.5, 
so we shall adopt the corresponding limits in this section. 



3.1 n 



This process occurs in the model via the one loop diagrams shown in Fig. 4. There are two 
couplings which are responsible for this process: and hi. In fact, the predictions of this 
model are similar to the ones discussed in Ref. [12], with one difference that here we have 
interference between diagrams generated by gij and those induced by hi. In the present 
model, ignoring the electron mass, which is an excellent approximation, the branching ratio 
is given by 



BR(/i -^erf) 



27a 
16ttG 2 f 



F(x d f^ + H(x u jV: 4 V a ^ 



(35) 
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Figure 4: One loop diagrams leading to \i — > e^y decay. 



where Xd t = m^./M| and x Ua = m^/Mf. The dimensionless functions -F(x) and H(x) are 
given by [26, 27] 



F{x) 



x (1 - x)(5 + x) + 2(2 x + 1) In x 
12 (1 - x) 4 ' 



„, . 1 (l-x)(5x + l)+2x(2 + x) lnx 

i] X = ; r-j . (OO) 

v ; 12 (1-x) 4 V ; 

The branching ratios for other £i — > processes can be derived in a similar way. The 
resulting constraints on the model parameters are summarized in Table 1. Here all of the 
experimental limits are taken from Ref. [20] except for /i — > limit which is taken from 
Ref. [25]. 

An interesting feature of this analysis is that the gry couplings are only weakly con- 
strained from these processes. This is owing to a GIM-like cancelation for the amplitude 
for this process from the first two diagrams of Fig. 4. This is similar to the model discussed 
in Ref. [12]. This cancelation occurs, in the limit of down quark mass being zero, since the 
charge of the internal leptoquark (2/3) is twice as large and opposite in sign compared to 
the charge of of the internal down quark (—1/3). The amplitude for the diagram when the 
photon is emitted from the scalar line is twice smaller compared to the diagram where it is 
emitted from the fermion line, which leads to the cancelation. The amplitude that survives 
has a suppression of (m^/M|g), which causes the weak limit. Because of this cancelation, 
we can derive correlated limits on the masses of the leptoquarks and the vector-like quark 
from jj, — y ey, since only the hi couplings are involved. This is shown as a contour plot 
in Fig. 5. To get the largest possible masses, we set the Yukawa couplings hi = hi — 1, 
as large as allowed by perturbativity. If /i — > e^y is discovered at the current limit [25] 
the masses should lie to the left of the red contour in Fig. 5, while a measurement of 
BR(,u — > ej) = 1.0 x 10 -12 would require the masses to lie to the left of the blue contour. 
The LHC reach for a leptoquark of this type is 1.5 TeV [28], which would serve as a cross 
check in this case. 
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Process BR limit 



Constraint 



fi ->• e 7 < 2.4 x 10" 12 
r -»■ e 7 < 3.3 x 10~ 8 
r -)■ ^7 < 4.4 x 



Table 1: Constraints on model parameters from ^ — >■ ^7. 

14 
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Figure 5: The allowed region for for the leptoquark mass Mi versus the vector-like quark 
mass M from fi — > ej. Here regions left of a contour is allowed for a fixed value of the 
branching ratio. Thus, if Br(fi — > ej) = 1.0 x 10~ 12 is measured, the masses must lie to 
the left of the blue contour. 

3.2 ti -> 3e 

In this process, the photon can be off-shell, and therefore, there is no GIM-like cancelation 
for the Qij contributions. It turns out that in addition to the photon penguin diagrams, there 
are also Z penguin diagrams and box diagrams (the Higgs boson exchange is suppressed 
by the electron mass). The mixing between vector-like quark and SM chiral quarks also 
plays a role in this process. The expression for the decay width is rather lengthy, which 
we do not present for brevity, but it is similar to the one given in Ref. [12]. To simplify 
the problem in deriving the constraints, we assume that only the top quark mixes with 
the vector-like quark, or equivalently s^, S24 ^ ^ 3 i n Eq. (14), where A ~ 0.22 is the 
Wolfenstein parameter, while S34 could be as large as 0.3, consistent with constraint from 
Z — )■ bb constraint [29]. For to 2 / 3 exchange (corresponding to down- type quark inside the 
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F(x di )^+H(x Ua )V: 4 V a4 ^ 



< 



6.6x10- 
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Figure 6: The tree level FCNC diagram leading to D — D mixing. 

loop), we assume that there is no accidental cancelation among the different couplings gij, 
and thus omit terms such as g\ 3 g 2 zgjk with j, k = 1,2. For degenerate leptoquark masses 
of 400 GeV, and for the vector-like quark mass set equal to 600 GeV, we obtain: 

\hh 2 \ < 2.7 x 10- 4 (3.4 x 10- 4 ); \g 13 g 23 \ < 1.7 x HT 3 ; \ 9lj g 2j \ < 8.6 x 10" 4 , j = 1,2, (37) 

for C34 = 0.98(1.0). These limits are obtained by assuming that contributions from one 
type of coupling dominates at a time. While these limits are stringent, they do not pose 
any restriction on the neutrino masses and mixings. The decay [l — > 3e may be within 
reach of next generation experiments, with the couplings lying in the range (10~ 2 — 1) and 
the leptoquark mass around a TeV. 

3.3 fi, — e conversion in nuclei 

Since this model features direct interactions of quark and lepton via the leptoquarks, fi — e 
conversion in nuclei occurs. The diagrams are similar to the ones discussed in Ref. [12], with 
tree level and loop contributions. There is a more direct link between neutrino mass and 
the loop induce fi — e conversion process. If we assume that only the top quark mixes with 
the vector-like quark as in the case of /i — > 3e, then there is no tree level a; _1//3 exchange 
contribution to this process. Following the procedure outlined in Ref. [12], from the limit 
on /i — e conversion in 48 Ti, we obtain (for Mlq = 400 GeV, and a vector-like quark mass 
of 600 GeV) 

\h x h 2 \ < 2.2 x 10~ 4 (9.8 x 10~ 3 ); l^s^sl < 8.7 x 10~ 4 ; \g n g 21 \ < 4.6 x 10" 6 , (38) 

for C34 = 0.95(1.0). Again, this analysis suggests that for natural values of the model 
parameters, this process may be within reach of next generation experiments. 

3.4 Tree level D° - D° mixing 

The FCNC that occurs in the up-quark sector (see Eq. (12)) induces tree level D° — D° 
mixing mediated by the Z boson, as shown in Fig. 6. The neutral Higgs boson induced 
contribution from Eq. (13) is suppressed by light quark mass and can be ignored. The 
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mixing amplitude is given by [30] 



Ml 



12 



new 



V2G F 
3 



(VuV 24 ) 2 m D f D B D r) D (n), 



(39) 



where m D = 1.9 GeV is the -D meson mass, /d = 0.201 GeV is the D meson decay constant, 
Bd(h ~ ttid) = 0.865 is the bag parameter, and t]d{h ~ mo) = 0.78 is the QCD correction 
factor. All the numbers here are taken from Ref. [31]. By using Am D = 1.6 x 10~ 14 
GeV [20] , one obtains the constraint 



According to Eq. (14) this constraint implies IC14S14S24I < 2.5 x 10~ 4 . As a result of 
this limit, unlike in a four generation model where there is no such FCNC process, the 
vector-quark contributions to meson mixing (eg., in the Bd sector) cannot be too large. 

3.5 B s — > fi + fi~ decay 

Recently, the CDF collaboration has reported a hint for the decay B s — > with 
BR = (l.8^o'g) x 10~ 8 [32]. LHCb collaboration has not confirmed such a hint, and quotes 
an upper limit BR(B S — > /i + /i~) < 1.4 x 10~ 8 [33]. The SM prediction for this branching 
ratio is BR(B S — > /i + /i~) = 3.2 x 10~ 9 , which means there is ample room for new physics 
in this process. 

Tree level exchange of leptoquarks can contribute to B s — > decay in our model. 

The branching fraction is given by [34] 



where the SM contributions have been ignored. This enables us to fit BR(_B S — > ) = 
1.8 x 10~ 8 , with M 2 = 400 GeV for the leptoquark mass, and | #22023! ~ 4 - 2 x 10 ~ 3 - 

3.6 New CP violation in B s — B s mixing 

Our model of leptoquarks and vector-like quark generates new contributions to B s — B s 
mixing. There are two sources, one through LQ induced box diagrams, and the other 
through SM-like box diagrams with the vector-like quark (see Fig. 7). Including these 
contributions, the B s — B s mixing amplitude becomes 



M 12s = [(V; 2 V 33 ) 2 V 33So(x 3 ) + 2 (V* 2 V 33 ) (V: 2 V 43 ) V 34So(x 3 , x 4 ) + (V 4 * V 43 ) 2 V uS (x 4 )} 



V U V 2 * 4 \ < 2.5 x 10 



-4 



(40) 




(41) 




(42) 
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Figure 7: Box diagrams leading to B s — B s mixing 



The functions So(x a ), So(x a , x@) with x a = vr? u /m^ are the Inami-Lim functions [35], 
whereas %-,^b with i,j = 3,4 are the QCD correction factors. The numerical values for 
these factors for a 600 GeV for vector-like quark mass are [36] 



f?33 = Vb = 0.5765, r] 34 = r] M = 0.514. 
It is sometimes more convenient to parametrize M 12s as [37] 



M 



12* 



Mlf s 



Mtf s 



r ls e i2ai ° + r 2s e t2(72 ° 



(43) 



(44) 



where {r ls ,o"i s } and {r 2s ,cx2 S } are the new contributions. With this parametrization, one 
can write 



AfTif( s 



Am 



SM 



1 + r ls e l2cTls + r 2s e 



s 



j/il>4> 



sin [2/3 s SM - Arg (l + r ls e l2(T1 ° + r 2s e i2a2 °)] , 



(45) 



with /3f M = Arg [(-V^V^,) / (^22^2*3)] = °- 019 ± °- 001 > Am B? = ( 19 - 3 ± 6 - 74 ) PS -1 - 

The main reason to highlight this phenomenon is because there are hints for new sources 
of CP violation beyond the SM in the D0 measurement inferred from the charge asymmetry 
in the same sign di-muon decay of the B mesons [16]: 



A, 



AH 



-(0.787 ±0.172 ±0.093)%. 



(46) 



N++ + jy 

Here N ++ (N ) is the numbers of events containing two b hadrons that decay semilepton- 
ically into two positive (negative) muons. Eq. (46) can be written as a linear combination 
of two asymmetries [16, 38] 



A b sl = (0.506 ± 0.043)4 ± (0.494 ± 0.043)a 



s 



(47) 
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where a 9 sl (q = d, s) is defined as [16] 



a 



<? _ 

si — 



T(B q -» fi+X) - T(B q -» /i-X) 
r(fi 9 ->• /i+X) + T(B q -> //-X) 



(48) 



In the SM, of, = -4.8l};!j x 10" 4 and a s sl = (2.1 ± 0.6) x 10~ 5 [39], so that {A b sl ) SM = 
— 2.31q6 x 10~ 4 which is 3.9a away from the current measurement (see Eq. (46)). A likely 
explanation is that there is a new source of CP violation in B s — B s mixing. 

Additionally, the measurements of relative phase between B s mixing amplitude and 
B s — > J/ip(f) decay amplitude (Sj/^ = sin 2/?/ // ^) as well as the measurements of decay 
width difference performed by CDF [40] and D0 [41] yield [42]: 



Here there is a 2.1<r discrepancy from SM prediction for /3 S , which may be another hint for 
physics beyond the SM. 

It is interesting to see whether the vector-like charge 2/3 quark can resolve these prob- 
lems. The best fit to the data for a fourth generation quarks, including the preferred values 
of the CKM mixing angles, is given in Ref. [43] which shows f3 s = 0.03 corresponding to 
ri s = 0.02 which is still far from the experimental central value (see Eq. (49)). This result 
should hold in the present model as well. We conclude that the mixing with vector-like 
quark is not enough to get the central value of (3 S , so the LQ induced box diagram is a more 
promising source for the new physics here. 

Ignoring the effect of extra family mixing, the LQ contribution {r2 S ,a2 S }, can satisfy 
the best fit given in Ref. [37], i.e. {0.5,120°}. This would correspond to Ig^tol ~ 0.05 
where the index % is summed. However, since the limits BK(B S — > fJ> + fi>~) ~ 10~ 8 and 
BR(5 S — > e + e~) < 2.8 x 10~ 7 must be satisfied, the dominant contribution should arise 
with the r lepton inside the box diagram loop. The phase is not constrained and therefore 
can take the fitted value of 120°. 

3.7 sin 2/3 versus ex 

Analogous to the B s system, this model also provides new contributions to B d — B d mixing, 
arising through LQ induced box diagrams and new mixing between SM quarks and the 
vector-like quark. Owing to the stringent constraint arising from the tree level D° — D° 
mixing, see Eq. (40), analogous mixing in the K° system would be more suppressed. 
The new contributions to the Bd mixing may explain the apparent discrepancy between 
sin 2/3 and determinations within the SM. The fitted value of sin 2/3 obtained from |e^|, 
Arrid/ Am s , \V u b/V c b\, BR(_B — > tu), and 7 measurements is given by 



pjm = o.47+°;^ u 1.091°; 



(49) 



sin 2(3 



.fit 



0.79 ±0.039 



(50) 
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Figure 8: The unitarity triangle fit leading to sin 2/5 determination. 

which differs by 3.1a from the world average experimental value [42], 

S J/i!Ks = sin 2/3 cxp = 0.671 ± 0.023. (51) 

The three family SM fit to the various CKM observables is shown in Fig. 8. A possible 
explanation is that there is new physics that affects the B d system, which we parametrize 

as 

M — M SM 

M su = ru^ + t^, (52) 

M 12d 



where 
M 12d 



Glm 2 



F" l W 



(V 3 \V 33 ) 2 m3So(x 3 ) + 2 (V 3 \V 33 ) (V&V43) %4^o(x 3 , x 4 ) + (V£V 43 ) 2 m S (x^) 



12vr 2 

+ Msf) (53) 

This is analogous to the discussion of B s mixing. With this formula, one can write 

S JN>Kt = sin [2/3 fit + <P B *] , (54) 

where 

<P B * = Arg [1 + r ld e i2aid + r 2d e i2a2d ] . (55) 

Unlike in the B s system, the effect of the vector-like quark can induce a significant effect 
to resolve the tension in sin 2/3 determination. In order to see this effect, let us ignore for 
the moment the LQ contributions. Then, if we choose [43] 

V 3 \V 33 = 0.009e ia56 

V£V4s = 0.00096e" a - 35 , (56) 
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Figure 9: Diagrams leading to the neutrinoless double beta decay through vector-scalar 
exchange. 

we obtain Sj/^k s = 0.68, which is in the good agreement with the experimental value. 
Note, that the LQ contribution is not strongly constrained by neutrino mass nor lepton 
flavor violation, so it can get close to the experimental value, as long as the LQ mass is less 
than 1 TeV. 



3.8 Neutrinoless double beta decay 

Although this model can accommodate inverted hierarchy in which neutrinoless double 
beta decay (/3/3q v ) may occur with sizable effect, it is still interesting to see that even in 
the normal hierarchy case, such process may be observed through vector-scalar exchange 
mechanism [18], depicted in Fig. 9. 

The effective Lagrangian of the new v — e c — u — d vertex after Fierz rearrangement is 



y*»ncw 



G p 
4\?2 



where 



u (1 + 75) d v e (1 + 75) e c 



+ -jua^ u (1 + 75) d VeP^ 



;i + 7 5 )e c 



Ml 



(57) 



(5f 



2V2MfG F 

This process is similar to MSSM models without i?-parity violation discussed in [18] and [44]. 
Following Ref. [44], and by using the results from Heidelberg- Moscow experiment on /3/3o^ 
decay rate [45], one obtains for Mi = 300 GeV and M 2 = 350 GeV 



\g!ihiVu\ < 4.3 x 10 



-7 



(59) 



The mixing matrix elements | V14V24 1 is constrained by D — D mixing process and has to 
be less than 10~ 4 . Since the coupling gn is not constrained by neutrino mass, it could be 
of order one. For Vu 10~ 5 and hi ~ 10~ 2 (assuming NH case and vector mass of order 
sub- TeV) from lepton flavor violation constraints, one ses that neutrinoless double beta 
decay might be observable even in the case of normal mass hierarchy. Of course, for this to 
be valid, the leptoquarks and vector-like quark have to be light. 
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4 Conclusions 

In this paper we have presented a new two-loop neutrino mass generation model which has 
the effective operator 3 of Eq. (2). Generating this effective operator in a renormalizable 
theory would require the addition of a charge 2/3 vector-like quark and a scalar leptoquark 
doublet tot the standard model spectrum. We have studied the phenomenology of this 
model. This model can explain the CP violation parameters in the B s and the Bd system. 
The leptoquarks of the model generate new CP violating contributions in B s — B s mixing, 
which can explain the di-muon anomaly reported by D0The apparent tension in the de- 
termination of sin 2/3 from B d decays and from the global analysis including from the 
K meson system also finds a natural explanation in this model. Neutrinoless double beta 
decay may occur through vector-scalar exchange and may be observable even with a normal 
hierarchy in the neutrino masses. 
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